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We investigate the effect of coupling (intensity and nature), applied field, and anisotropy 
on the spin dynamics of a multi-layer system composed of a hard magnetic slab coupled to 
a soft magnetic slab through a nonmagnetic spacer. The soft slab is modeled as a stack 
of several atomic layers while the hard layer, of a different material, is either considered 
as a pinned macroscopic magnetic moment or as an atomic multi-layer system. We com- 
pute the magnetization profile and hysteresis loop of the multi-layer system by solving the 
Landau-Lifshitz equations for the net magnetic moment of each (atomic) layer. We study the 
competition between the intra-layer anisotropy and exchange interaction, applied magnetic 
field, and the inter-slab exchange, dipolar or Dzyaloshinski-Moriya interaction. Comparing 
the effects on the magnetization profile of the three couplings shows that despite the strong 
effect of the exchange coupling, the dipolar and Dzyaloshinski-Moriya interactions induce a 
slight (but non negligible) deviation in either the polar or azimuthal direction thus providing 
more degrees of freedom for adjusting the spin configuration in the multi-layer system. 

PACS numbers: 75.10.-b General theory and models of magnetic ordering - 75.78.-n Magnetization 
dynamics - 75.10.Hk Classical spin models 



I. INTRODUCTION 

Although the "nano-rush" tends to dominate the realm of technological applications, especially 
magnetic recording, multi-layer magnetic systems benefit from a growing interest in this area mainly 
due to their high performance |THH]. On the other hand, magnetic multi-layer systems and thin films 
benefit from the acquired long-standing experience and know-how both in growth and characteriza- 
tion leading to a good control of the relevant intrinsic parameters (dimensionality and anisotropy). 
Moreover, there are many well-established techniques for precise measurements, such as FMR [1HS], 
BLS [10], and the ever developing optical techniques |11H14| . to cite a few. 

The magnetization dynamics of laterally confined elements of alternating magnetic and nonmag- 



* Electronic address: 



andres.franco@univ-perp.fr, hamid.kachkachi@univ-perp.fr 



2 



netic layers exhibits a large variety of interesting phenomena for both applications and fundamental 
research. In this context, theory has to play its usual role of providing reasonable models for inter- 
preting the observed phenomena and suggesting new experiments. An issue of particular interest in 
this context concerns the coupling in multi-layer systems [15J. The inter-layer coupling determines 
the mechanisms of transport and propagation of a stimulus applied at one end of the structure and 
the mechanism of the adjustment of magnetic configurations in metallic multi-layer systems by an 
external magnetic field. In the context of perpendicular magnetic recording with exchange-spring 
media, improved write-ability is achieved by appropriately tuning the coupling between the soft 
and hard magnetic layers |2"| |16I|2"U] , Therefore, it is of paramount importance to fathom the nature 
of interaction acting at the interface and the role it plays in conveying any perturbation through 
the multi-layer system. There is a great amount of published work investigating the inter-layer 
coupling and its effects on the magnetic properties [8j |21-28j. For instance, in Ref. El the authors 
investigated the magnetization dynamics due to spin currents in magnetic double layers and argued 
that transport in this structure is governed by a kind of long-range interaction called the dynamic 
exchange coupling. The microscopic mechanism underlying this effective coupling and the way it 
affects the collective behavior of magnetic hybrid structures require further investigation. While 
the RKKY coupling provides an interpretation of experiments on very thin conducting spacers as- 
suming a high degree of perfection [29], rough surfaces may induce strong stray fields and thereby 
dipolar coupling in multi-layer systems. For a magnetic bi-layer several configurations are consid- 
ered in Refs. [22j [3011321 In particular, in Ref. |22] the calculations confirm the fact that in the 
absence of roughness the dipolar coupling between two perfectly flat infinite planes vanishes and 
that for planes of finite dimensions the dipolar coupling may give rise to a ferromagnetic or an 
anti-ferromagnetic coupling. On the other hand, one cannot exclude the Dzyaloshinski-Moriya in- 
teraction, which is an anti-symmetrical exchange interaction and mainly stems from a combination 
of low symmetry and spin-orbit coupling |33| I34j . In the presence of disorder, especially at the 
interface of thin films or multi-layer systems, the Dzyaloshinski-Moriya interaction has been shown 
to play an important role since local symmetry is broken by surface effects. Indeed, it leads to 
large anisotropy and may even change the magnetic order, see Ref. f33B36] . In particular, it has 
been shown that the Dzyaloshinski-Moriya interaction is induced by spin-orbit coupling between 
two ferromagnetic layers separated by a paramagnetic layer |23] , 

In the present work, we bring a new contribution that attempts to further clarify the role of three 
inter-slab couplings, namely exchange, dipolar and Dzyaloshinski-Moriya, in the spin dynamics 
and magnetization profile (MP). It is also essential to compare these interactions with respect 
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to their efficiency in the realignment of the spin configurations and eventually the magnetization 
reversal. More precisely, we investigate the effect of coupling (intensity and nature), applied field, 
and anisotropy on the spin dynamics of a multi-layer system composed of a hard magnetic slab 
coupled to a soft magnetic slab through a nonmagnetic spacer. The soft slab is modeled as an array 
of several atomic layers while the hard layer, of a different material, is either considered as a pinned 
macroscopic magnetic moment or as an array of atomic planes. We compute the magnetization 
profile and hysteresis loop of the whole system by solving the (coupled) Landau-Lifshitz equations 
for the net magnetic moment of each (atomic) layer. 

This paper is organized as follows. After the introduction we define the system we study in 
section II. In section III we deal analytically and numerically with the rigid interface. Section 
IV deals with the soft interface. It presents the results for the effects of applied field, inter-slab 
couplings, and their comparison. It ends with the results of hysteresis loops. We finally summarize 
the main results in the conclusion and give a few perspectives for future investigations. 



II. STATEMENT OF THE PROBLEM 



We consider the coupling between two slabs: one is a sufficiently thin hard magnetic system 
(HMS) with strong out-of-plane anisotropy and, on top of it, a thick soft magnetic system (SMS) 
with in-plane anisotropy. The SMS slab will be modeled as a multi-layer system of atomic planes of 
e.g. Fe, whereas the HMS slab will be modeled either as i) a single macroscopic magnetic moment 



in which case we have a rigid interface (RI), see Fig. 1(a) or as ii) a multi-layer system of atomic 



planes of e.g. FePt, and in this case we get a relaxed or soft interface (SI), see Fig. 1(b) In both 
cases, i.e. SMS /HMS-RI or SMS/HMS-SI we have an exchange-spring system with a characteristic 
MP that depends on various physical parameters such as the SMS-HMS coupling, the anisotropy 
and exchange coupling within each slab, and the applied field. Accordingly, in this work, we will 
investigate the effects of these parameters on the MP in the whole system and on the hysteresis 
loop for the two configurations. In particular, we will consider and compare the effects of three 
kinds of inter-slab couplings, which are exchange, dipole-dipole and Dzyaloshinski-Moriya. For 
later reference, we denote by EI-ES, DDI-ES, and DMI-ES the exchange-spring system where SMS 
and HMS are respectively coupled by exchange interaction (EI), dipole-dipole interaction (DDI), 
or Dzyaloshinski-Moriya interaction (DMI). 

In Ref. |37] the relatively simpler case was considered where both the SMS and HMS were rep- 
resented by their respective net magnetic moments. The relaxation processes of the corresponding 
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(a) Rigid interface 



(b) Soft interface 



Figure 1: Scheme of a hard/soft coupled bi-layer system 



magnetic dimer were then studied. In particular, the relaxation rates related with the reversal of 
its magnetization were computed. Moreover, the effects of the three couplings mentioned above 
were investigated for different configurations of the anisotropy easy axes, and at finite temperature. 
In the present work, we are interested in the exchange-spring effect and in the MP at equilibrium 
and at zero temperature. Therefore, the main objective here is to investigate how the magnetic 
structure/state of the SMS adapts to a change in the physical parameters, especially at the interface. 



III. RIGID INTERFACE : NUMERICAL VERSUS ANALYTICAL CALCULATIONS 
A. Magnetization profile of the multi-layer system 

The exact solution for the MP has been obtained for the system with RI, in the absence of the 
external magnetic field |38j . For definiteness, the SMS-HMS lies on the xy plane and the out-of- 
plane anisotropy in the HMS slab points in the z direction. We denote the (net) magnetic moments 
of HMS by <r, and the M s magnetic moments of the SMS by Sj,i = l...,M s . Therefore, there is a 
total of M s + 1 magnetic moments in the system pictured here as an open chain of coupled magnetic 
moments. Indeed, it is assumed that all atomic magnetic moments of each sub-layer i in the SMS 
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are parallel to each other so that their magnetic state can be represented by the net (normalized) 
magnetic moment Sj. Hence, we do not include boundary effects in this study. For the HMS, cr 
represents the net magnetic moment of the whole slab. The system setup is shown in Fig. [2] 
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Figure 2: Setup of the SMS/HMS system with rigid interface. 

If we denote by the anisotropy easy axis (along the z axis) of HMS and e s the anisotropy 
easy axis, taken here along the x axis, of the SMS, the total energy of the system can be written as 

E = -EsMS + ^HMS + -^Int, 

with 

E S ms = -D s ( s f ) 2 - J s X) Si ' Si+1 

j=l i=l 

being the energy of the SMS slab. The first term is the in-plane anisotropy with constant D s (s 
stands for soft), the second term is intra-SMS nearest-neighbor exchange coupling, of intensity J s . 
Similarly, the energy of the HMS slab contains the out-of-plane anisotropy contribution and the 
exchange coupling. However, in the present case of a RI, the energy coupling is dropped since we 
assume that all atomic magnetic moments within the slab are already tightly bound together by 
the exchange coupling. On the other hand, the anisotropy contribution is also dropped because it 
is assumed to be strong enough to pin the macroscopic magnetic moment cr in the z direction. 

Finally, the interaction between the two slabs will be denoted in the sequel by A and will be of 
three different origins : exchange, dipolar or Dzyaloshinski-Moriya. In the present case, A is the 
exchange coupling and will be set to Jq, so that 



6 



Eint = -Act • si = -Jo<r • si. 

The analytical solution of the problem requires an additional simplification that consists in 
setting Jo = J s = J. Therefore, we can write the energy (in units of exchange coupling) as 

i=l i=l 

where we have introduced the dimensionless parameter d s = D s /J s . 

Then, we introduce the coordinate z of each atomic sub-layer along the z direction, such that 
z = and z = L represent the sub-layers at the interface and the top of the SMS slab, with 
L being the thickness of the SMS slab and a that of its atomic sub-layers, i.e. L = a(J\f s — 1). 
Note that i = represents the HMS slab and a slight shift of 0.5a in the reference frame between 
the numerical and the analytical calculations should be taken into account. Indeed, the analytical 
calculations were performed in the continuum limit. Moreover, the angle deviation with respect 
to the z axis 8i — 6i+\ of two adjacent magnetic moments Sj,Sj+i is supposed to be small. For 
convenience, we denote £j = 9i — 9 e , where 8 e is the equilibrium polar angle of the multilayer system 
as defined in [38J. For typical materials d s <C 1, thus 9 e ~ 0. In such a case one can write down the 
(finite-difference) equation for £ from the energy minimization. The solution of this equation can 
be written as |38j 

r£ drj 

Z ~ JO VC L +ds cos 2jj' 

(2) 

C L = (d s sin2££,) 2 - d s cos2£c, 

£l being the angle deviation of the SMS top sub-layer, i.e. for z = L and is given through the 
relation 

L=t- d l . (3) 

J VC L + d s cos2£ v ; 

o 

There is a minimal number of sub-layers of the SMS slab, A/" s mm , or length of the chain, L min , 
necessary for an onset of noncolinearities of the magnetic moments Sj. In Ref. 38, 7V™ m was found 
to be given by 

A4 min ~ — J=. (4) 



This yields, for instance, A/" s mm ~ 11 for ds = 0.01. 
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B. Magnetization profile of the exchange spring: numerical solution 



In order to obtain the MP for the system described above, we solve the set of coupled Landau- 
Lifshitz equations (LLE) 



1 dsi 

7~dt 



s,- x h 



cff 



as, x ( s,- x h! 



off 



(5) 



where hf ff 



-Se/Ssi is the (dimensionless) effective field comprising the anisotropy and exchange 
contributions, together with the interaction at the interface, a is the phenomenological damping 
parameter, set here to 0.01. In these zero-temperature calculations, damping is used to drive the 
system into the equilibrium state, i.e. the state of minimal energy. We start from a state of 
homogeneous magnetization along the HMS anisotropy axis and allow the system to relax towards 
a minimum that yields the MP of the system. 

In Fig. [3] we plot the angle deviation £J of the magnetic moment of the top SMS sub-layer 
against the number of sub-layers in the SMS slab. The results for three values of the in-plane 
anisotropy d s show a good agreement between our numerical results and Eq. (|3j). 
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Figure 3: Comparison between analytical and numerical calculations of the magnetization of the top sub- 
layer £^ for various values of the SMS in-plane anisotropy, d s = 0.01 (left), d s — 0.025 (middle), d s = 0.05 
(right). 

One can also see how J\f m - m depends on d s through Eq. Q . For a weak d s the system behavior 
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is dominated by the fixed magnetic moment of the HMS slab through the interface exchange cou- 
pling. Adding more sub-layers to the SMS slab attenuates this effect since the interface interaction 
is counter-balanced by the anisotropy and intra-layer interaction of each SMS sub-layer, leading 
eventually to a deviation of the magnetic moment of the top SMS sub-layer for J\f > A/" m ; n . On the 
other hand, when d s increases, the contribution of each SMS sub-layer to the energy of the whole 
system increases, thus leading to a smaller A^min- 

IV. SOFT INTERFACE 

In the previous section, we defined a model for a system with rigid interface and energy given 
by Eq. ([!]). The agreement with the previously published analytical work validates the numerical 
approach used here, which can now be employed to deal with more complex systems that are 
difficult to tackle analytically. 

More precisely, we extend the previous model by replacing the macroscopic magnetic moment 
<j of the HMS, that was pinned along the anisotropy direction, by Ah sub-layers each with an 
out-of-plane (in the z direction) uniaxial anisotropy of intensity D^. Moreover, in the present case 
we apply an external magnetic field H at an angle Oh with respect to the z axis. Now, within the 
HMS slab, and similarly to the SMS slab, each sub-layer of index k carries a magnetic moment er^ 
with k = 1, . . . ,J\fh, see Fig. [|] 

The energy of this model then reads 

E = -EsMS + ^HMS + -^Int, 

with 

Esms = ~VsH ■ e H - D s ^2 ( s f f ~ Js ^2 Si- s i+1 

i=0 i=l i=l 

which is the energy ([I]) to which we have added the Zeeman contribution. 

E-rms = —fJ-hH <T fc • en — Dh ( a k) 2 ~ Jh ""fe-Cfc+i- 

fe=l k=l k=l 

Us and /i/j are the atomic magnetic moments of the soft and hard materials, respectively. 

The interaction contribution is a function of si and crj\f h and its form depends on the nature 
of the coupling between the two slabs, i.e. the top layer of the HMS slab and the bottom layer 
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Figure 4: Setup of the SMS/HMS system with soft interface. 

of the SMS slab. Here again, we use the same normalization as in Eq. ([I]) and divide by the 
exchange coupling J s within the SMS slab, and introduce the following dimensionless parameters 
(with obvious notation). 

Hence, for the interaction energy we write 



£int = AJ 7 (si,o- A ^) 



where for exchange we have 



for DMI 



and for DDI 



£int = — Aei^a/; • si, Aei 



Jhs 

X' 



£lnt = — ^DMI e DMI " ( a Mh x s l) J ^DMI = 



D 



4vr J, 



£lnt = ^DDlCNh ■ ^12Sl, AdDI 



(6) 



(7) 



(8) 



with 



Pi2 = 3(^ .e a e a - -►)-!. 



(9) 
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being the DDI tensor. 

h = fj, s H/J s is the reduced magnetic field; Jh s = Jh/Js is the inter-layer exchange interaction 
between the HMS and SMS, D is the magnitude of the DM vector, and d is the distance between 
the SMS and HMS. Furthermore, we assume that the magnetic moment of each sub-layer, within 
each slab, has the same magnitude. 

Now, we are ready to discuss the results of several numerical studies investigating the effects 
of the applied magnetic field h, the in-plane anisotropy d s of the SMS slab, and the nature of the 
inter-slab coupling and its intensity. In each case we compute the magnetization profile, namely 
the angle deviation with respect to the vertical (z) axis of each layer as we move upwards from the 
bottom layer of the HMS slab to the top layer of the SMS slab. We also compute in each situation 
the angle deviation of the bottom layer of the HMS £l and the top layer of the SMS slab (g, as a 
function of the applied field. 

The parameters used for the calculations are dh = 0.02, = 1.44 and Mh = 10. Note that 
the change from HMS to SMS occurs when the sub-layer index n passes from the top HMS layer 
to the bottom SMS layer, with n = 0,...,Mh + A4 — 1 being the sub-layer index that labels the 



sub-layers of the whole system (SMS+HMS). M s = 41 for the calculations presented in section IV A 



and M s = 11 for those presented in sections |IV B 1 to IV B 3 In order to better understand the 
effects of each variable we set to zero all those that we are not studying, with only two exceptions, 
namely d s , which takes the value of 0.01 for all the profile plots and various values for (^ and (ft 
variation. Aei = 1.44 for calculations with varying h and 9h, but zero otherwise, except when we 
study the effect of a varying exchange inter-slab coupling. The values of the parameters considered 
here are typical of the Fe/FePt multi-layer systems [see Ref. [38land references therein]. 



A. Effects of the applied field and in-plane anisotropy 

The angle deviations, (ft for the HMS bottom layer and ( S L for the top layer of the SMS slab, 
are computed for various values of the applied field h with a direction 0h = and Aei = 1-44; the 
results are shown in Fig. [5] 

It is clear that the external field competes with the anisotropy and tends to align all the magnetic 
moments of the system parallel to its direction. Indeed, with enough layers in the SMS slab, (^ 
can be obtained from the Stoner-Wohlfarth equation 

h sin (6 -9 h ) -d s sin 26 = 0. (10) 
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Figure 5: (a) System setup with inter-slab exchange interaction and applied field along the z axis; (b) 
variation of the MP; (c) the deviation of the magnetic moment of the bottom HMS sub-layer (fi and (d) 
of the top SMS sub-layer £| with the magnitude of the applied magnetic field with 9h = 0. Afh = 10 and 
N s = 41. 



In the specific case of 9^ = 0, we get 9 = arccos (h/2d s ), in agreement with the asymptotes (dashed 



lines) in Fig. |5(b) Indeed, Figs. 5(c) and 5(d) show that there is a critical value of the field, i.e. 
h c = 2d s , at which all magnetic moments are aligned along the direction of the field, i.e. 9h = 0. In 
this the applied field is along the HMS anisotropy, h c is not affected by the latter. Hence, 

for a system with a sufficient number of layers in the SMS slab, h c will depend only on d s . In 
general, however, h c should depend on dh, d s , and the applied field orientation. In cases where the 
system has less layers than the necessary number to attain the asymptotic value of ^ and given 
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by the SW equation (10), the magnitude of h c should decrease. 
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Figure 6: (a) System setup with inter-slab exchange interaction and applied field of varying direction; (b) 
variation of the MP; (c) deviation of the magnetic moment of the bottom HMS sub-layer (ft and (d) of the 
top SMS sub-layer £| with the field direction 9 H . h = 0.01, TV,, = 10 and N s = 41. 



Fig. [6] shows the same plots for a variable field orientation and the SMS intra-plane anisotropy 
d s , for h = 0.01 and Aei = 1.44. We see that as the field is turned towards the SMS easy axis, 
the MP obviously shifts upwards with no noticeable change in shape, reaching a deviation at the 
top sub-layer that is again given by the numerical solution of the SW equation (Fig. |6(b)| - dashed 
lines). Figs. 6(c) and |6(d)"] suggest that beyond a given Oh the increase in the deviation is almost 
linear with a slope that depends on the anisotropy d s . This implies a constant rate of change of (ft 
as it is mainly affected by dh (constant here), whereas the rate of change of & varies with d s . 
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B. Effect of inter-slab coupling 

Now, we investigate the effect of the inter-slab coupling, considering successively EI, DDI, and 
DMI. In the end we compare their effects on the magnetization profile. 



1. Exchange interaction 




The results in Fig. 7(a) show that due to the strong anisotropy of the HMS slab, varying the 
inter-slab exchange coupling only affects the sub-layers in the SMS slab that starts here at n = 11. 
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On the other hand, the inter-slab exchange coupling competes with the intra-slab exchange coupling 
J s and anisotropy D s of the SMS slab. Thus the weaker Aei the stronger is the deviation. 



Figs. 7(b) and 7(c) show the deviation (^ and ( S L as functions of Aei, respectively. As it could 
be expected, with increasing EI we achieve higher £j and lower (^ deviations. As we effectively 
increase the rigidity of the interface, the HMS deviation at the interface thereby increases, whereas 
that of the SMS decreases. This change in deviation is then conveyed through all the sub-layers 
by means of the intra-layer EI, leading to the observed changes in (ft and ( S L . For low values of d s , 
the # 

increase is much slower than the ( S L decrease. However, as d s approaches dh = 0.02, the two 
become comparable as the two systems are close to being identical. 

2. Dipolar Interaction 

In the present work we consider the system setup where the dipolar coupling is assumed to 
induce a ferromagnetic coupling between the two magnetic slabs. 

In Fig. |8(b) are plotted the MP for different values of the DDI inter-slab coupling Addi- Apart 



from the obvious shift upwards as Addi decreases, there is an abrupt change at the interface espe- 
cially for small Addi, mainly due to the fact that the in-plane anisotropy d s has a stronger effect 
than DDI. 

Figs. 8(c) and 8(d)| show that the system tends towards an asymptote as Addi increases. In our 



case, since the DDI vector is parallel to the HMS anisotropy (along the z axis), a strong enough 
interaction aligns the magnetic moments at the interface (and all the magnetic moments in the 
HMS) in the direction 9 = 0, thus driving the system into an effective RI state. As such, the 
asymptote can be found by calculating using the analytical expressions 

[38], i-e. Eqs. (g and 

Q. Indeed, Fig. 8(e)| where the MP is plotted for different values of Af s with very strong DDI 



(Addi ~ 5), shows a perfect agreement. Furthermore, if we examine the analytic curve (black 
line) near n = 20 we observe a regime where a slight variation in M s induces a large change in 
indicating that J\f s = 10 is approximately the critical length of the SMS chain. A similar regime 
starts to be seen for stronger values of DDI in Fig. |8(b)| , Then, if we take into account that this 
behavior is not observed for weak DDI, it suggests that the critical length of the chain increases 
with increasing DDI. 

DDI is a long-ranged interaction that can penetrate through the SMS leading a priori to a 
coupling of the HMS interface with every SMS sub-layer, and vice-versa, with a strength ( ~ <i~ 3 
where d is the distance between two magnetic moments. We also show in Fig. [8] that with Addi 
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Figure 8: (a) System setup for a system with inter-slab dipolar interaction with a dipolar vector parallel 
to the HMS anisotropy; (b) variation of MP; (c) deviation of the magnetic moment of the bottom HMS 
sub-layer £^ and (d) of the top SMS sub-layer £| with the inter-slab DDI Addi- (e) Comparison between 
PJ analytical expressions and strong DDI numerical calculations. Mh = 10 and Af s = 11. 
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Figure 9: Comparison between different MP for a RI system with DDI at the interface only and with an 
effective value for both A D di and j s . X DDI = 0.15, A£f DI = 0.201, jf = 1.095. 

at the interface with the bond ecm along the z axis, we can approximate the configuration of 
our system by a RI. Upon taking long distance interaction into account, the HMS interface will 
be coupled with every SMS sub-layer, but the SMS interface will only be coupled with the HMS 
at the interface, due to the latter being in a RI configuration. We call this a long distance (LD) 
configuration. Fig.[9]shows a comparison between the MP for Addi = 0.15 only at the interface and 
LD configuration. The effect of the coupling penetration is a global decrease in the deviation of the 
magnetic moments of the SMS sub-layers. A similar behavior can be obtained with the effective 
couplings A^pj = 0.201 and = 1.095 with interaction only at the interface. This means that the 
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effect of LD configuration is equivalent to that of an interaction which is limited to the interface, 
but with the re-normalized couplings A^j and Figs. [9] (a)-(c) show that A^-,j and do not 
change with the number of SMS sub-layers. However, Fig. [T0| shows that they do depend on Addi- 




Sublayer index Sublayer index 



(a) Addi = 0.1, Addi = 0.1285, jf = 1.0778 (b) Addi = 0.2, A D di = 0.2815, jf = 1.104 

Figure 10: Comparison between MP for a PJ system with DDI at the interface only, LD regime, and for the 
case with an effective values for Addi and j s yielding the same MP as for the LD regime. Af s = 10. 



3. DM Interaction 

As discussed in the introduction, it is relevant in the present study to investigate the effect 
of DMI on the dynamics of the MD, on the same footing as the (symmetrical) effective exchange 
coupling and anisotropic dipolar coupling. DM interaction may be induced by spin-orbit coupling 
between two ferromagnetic layers separated by a paramagnetic layer |23} 135] . It plays an important 
role in the presence of roughness and disorder at the interface of multi-layer systems. For a simple 
cubic lattice, on the (10 0) surface the DMI vector D lies in the layer plane and thus induces a 
perpendicular anisotropy. In the present work, we consider two situations where the DMI vector 
D lies either along the x direction (thus in the xy plane) or along the z direction {i.e. normal to 
the xy plane). 

a. DM vector in the x Direction 
According to Eq. ^fj, the DMI tends to orientate the magnetic moments in such a way that 



they are normal to each other and perpendicular to the vector D. Fig. 11(b) shows that with D 



along the x axis, the HMS magnetic moments at the interface align along their easy axis e z , as it is 
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perpendicular to the vector D and there are no other fields that could lead to a deviation from it. 
Therefore, all magnetic moments of the HMS slab will stick up along the direction 9 = 0. On the 
other hand, the SMS magnetic moment at the interface, because it has to be perpendicular both to 
the HMS magnetic moment {z axis) and the vector D [x axis), lies in the xy plane. However, in the 
SMS slab there is a competition between the DMI and the anisotropy d s in the xy plane leading to 
a variable azimuthal angle tp s as can be seen in Fig. [II] (c). With weaker DMI the anisotropy has a 
stronger effect and the SMS magnetization at the interface lies near the x axis. As DMI increases 
the magnetization turns towards the y axis as it has now to be perpendicular to the vector D. 



In Fig. 11(d) we see that the HMS magnetization stays along the z axis (£j = 0) while that 
of the SMS slab remains in the xy plane = 7r/2) for values of d s = 0.005,0.01,0.015 and 0.02. 
This means that the SMS magnetic moment at the interface is not affected by a varying SMS 



anisotropy in this setup. However, as is seen in Fig. 11(e) the (p 8 dependence on the SMS-HMS 
coupling is affected by the anisotropy, exhibiting higher deviations for lower values of the anisotropy. 
Indeed, it is more favorable for the interaction to win the competition and drive the SMS interface 
magnetization onto the y axis. For a strong enough interaction the SMS magnetization at the 
interface is pinned along the y axis and the system can be modeled as a rigid magnet where the 
magnetization of the SMS varies in-plane with no out-of-plane component. In this situation, the 
top SMS magnetization can be calculated using Eq. (|3j) with now the angle variation £l referring 
to ip s L instead, because £ s is fixed at ir/2. 



Fig. 12 shows the MP similar to the one shown in Figs. 11 (b) and (c), but in this case we 



force the HMS magnetic moments to have a deviation £ h = 7r, opposite to that of Fig. 11 (b), and 
we let the system evolve. This causes the orientation of the SMS magnetic moments to switch as 
well, thus demonstrating that for a system with this specific setup, it can be used as a "magnetic 
switch" since switching only the magnetic moment of the HMS slab induces a switching of the SMS 
magnetic moment. This is also true the other way round as is the case in exchange-spring systems 
where one attempts to achieve the reversal of the hard layer by smaller DC fields upon acting on 
the soft layer. 

b. DM vector in the z Direction 
Similarly to the previous case, the DMI will align the magnetic moment of one of the layers at 
the interface along the anisotropy axis of the corresponding slab, whereby its magnetization will 
not vary with the interaction. In the present case of D \\ e z , however, it is the SMS magnetic 
moment that remains constant because its anisotropy axis is normal to the vector D. On the other 
hand, the HMS (net) magnetic moment becomes pinned in the yz plane, implying that the angles 
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Figure 11: (a) System setup with inter-slab DMI interaction with a vector D perpendicular to the HMS easy 
axis; (b) variation of the polar MP (£); (c) azimuthal MP (if), (d) the polar and (e) azimuthal deviations of 
the magnetic moments at the bottom HMS sub-layer £^ and the top SMS sub-layer £ S L with the inter-slab 
DMI Addi, for various values of d s . 
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Figure 12: Magnetic "switch" behavior in the (a) out-of-plane (£) and (b) in-plane (</?) directions of the 
system with DMI at the interface and vector D along the x axis. A//, = 10 and Af s = 11. 

</? for both layers will remain constant and equal to tp h = ir/2 and tp s = 0. As such, the deviation at 
the interface for the HMS slab will vary according to the competition between the HMS anisotropy 
and the DMI, see Fig. 13(b)| Hence, changes in the SMS anisotropy will not affect the MP of the 
system. Indeed, calculations of the £^ and for various values of d s were performed and the same 



results, shown in Fig. 13(c) , were obtained for all of them, thus confirming that the system is not 



affected by varying the value of the SMS anisotropy. 

For larger values of DMI the system can be viewed as an "inverted" rigid magnet, where the 
SMS and the HMS (net) magnetic moments at the interface are pinned in-plane. The out-of-plane 
variation of the HMS magnetic moment results from the competition between the HMS intra-slab 
exchange interaction and anisotropy. In this situation, the variation in the HMS magnetization can 
be calculated by using Eq. (|3]), upon substituting dh for d s . 

4-- Comparison between the inter-slab couplings 



In Fig. [l4] we present different magnetic profiles for typical values of the EI (Xei = 1-44), DDI 
(Addi/Aej » 1(T 2 ) and DMI (Admi/Abj » 0.1), as can be found in Refs. [351 EH US] . Figs. [l4](left, 



right) present the MP in the polar (£) and azimuthal (cp) directions, respectively. In the polar MP 
the black curve with circles is the MP with only EI at the interface and serves as a reference. The 
red curve with squares is the MP obtained as we add DDI. Compared to the EI strength, the DDI 
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Figure 13: (a) System setup with inter-slab DMI interaction with a vector D parallel to the HMS anisotropy; 
(b) variation of the MP; (c) deviations of the magnetic moments at the bottom HMS sub-layer $ and the 
top SMS sub-layer £| with the inter-slab DMI A DD i. Af h = 10 and Af s = 11. 



is two orders of magnitude weaker and thereby its contribution to the alignment of the magnetic 
moments at the interface is very small. Nonetheless, it tends to align the interface along the z axis 
and this effect propagates throughout all sub-layers by the in-plane exchange interaction, causing a 
global decrease in the MP. The dark blue curve (triangles up) represents the MP for DDI only. We 
see that the weak DDI is not sufficient to overcome the anisotropy of each sub-layer, leading only to 
a slight deviation near the interface but away from the latter each slab remains parallel to its easy 
axis. We observe an increased (induced) deviation in the SMS as compared to that of the HMS. 
This is due to the stronger anisotropy of the HMS enhanced by the DDI-induced anisotropy in 
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the z direction, making the out-of-plane direction more favorable at the interface than the in-plane 
direction. On the other hand, for the reasons given earlier, the DMI with its vector D along the x 
axis, leads to a similar result with each slab aligned along its own easy axis even near the interface. 
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Figure 14: Effect of EI, DDI and DMI (x direction) on the magnetic profile of a exchange spring, with 
j = 1.44, Addi = 0.0144 and S = 0.144. Af h = 10 and JV S = 11. 



The MP in the azimuthal direction shown in Fig. 14 (right) compares the three inter-slab 



couplings. The EI and DDI are two interactions that do not induce an azimuthal rotation and the 
corresponding MP always lies in the same plane (92 = 0). The green curve (with diamonds) is the 
MP obtained after adding DMI with its vector D in the x direction. In Fig. [14] we see that the DMI 
induces a slight decrease in the deviation ip s because the DMI tends to align the magnetic moments 
of the SMS slab along the y axis, as is seen previously in Fig. ll(e)| Even though the magnitude of 
this interaction was not high enough to overcome the anisotropy and internal exchange, it induced a 
deviation at the interface in the azimuthal direction. This deviation is, however, constrained by the 
EI at the interface and the anisotropies of the slabs. The light blue curve (triangles down) represents 
the MP for DMI with no exchange at the interface. As it could be expected, the magnetic moments 
of the HMS are all aligned along the z axis (recall that D is along the x axis), while those of the 
SMS slab are all in the xy plane with a higher ip deviation at the interface than in the previous 
case. This is obviously due to the fact that now the DMI competes only with the anisotropy of the 
SMS slab. 
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5. Hysteresis loops 

Here we present a succinct study of the hysteresis loop for different values of the interaction. 
It helps us to understand how the switching mechanism of the multi-layer system changes with 
the nature and strength of the inter-slab coupling. All the hysteresis curves presented are plots 
of the normalized magnetization M z /M a versus reduced applied magnetic field h, along the HMS 
anisotropy axis (z axis). Mh = 12 and M s = 16. We start in a state of zero magnetization {i.e. the 
magnetic moments of all sub-layers are on the plane with no applied field) and iterate the LLE Q 
until we reach the state of minimal energy. Next, we increase the applied field slightly and wait 
for the system to reach the new equilibrium state. We keep increasing the applied field until we 
reach saturation of the magnetization. Then, we ramp the field down and up again thus closing the 
hysteresis loop. Obviously, for each value of h, we wait until the state of minimal energy is reached. 

Figure [l5] (a) shows hysteresis cycles for different values of EI. Let us denote by h c s the SMS 
switching field, by h c h the HMS switching field, and by /i sa t the saturation field of the whole multi- 
layer system. With no interaction (Aei = 0) and for h = the equilibrium state is that where the 
magnetic moments (of HMS and SMS) are oriented along their respective easy axes, i.e. 9k = 
and 9i = tt/2 for all k £ HMS,i G SMS. The net magnetization of the system is that given by the 
HMS net magnetic moment projected on the z axis. When h is increased the magnetic moments 
of the SMS start to rotate towards the applied field in a reversible process until the system reaches 
saturation {6k = 0i = for all k, i) at h = /i sat . As we ramp down the field across zero until it 
reaches the SMS anisotropy field {h c s = 2d s = 0.02), the SMS magnetic moments switch towards 
the field direction 9{ = ir for all i £ SMS. Further increase of the applied field (in the opposite 



direction) induces a slight reversible deviation of the HMS magnetization. In Fig. 15 (a) this 
corresponds to the plateau (on the Aei = curve) from h = —0.02 to h ~ —0.04, until h reaches 
the HMS anisotropy field hi = h Ba t — 2dh = 0.04. At this value of the field, the HMS magnetic 
moments coherently switch from 9^ ~ to Of. = tt, thus achieving negative saturation. Along the 
lower branch the system follows the same switching process from 6 = it to 9 = 0. 

For nonzero but weak coupling, Aei = 0.04, the system follows the same behavior, but neither 
the HMS nor the SMS goes through a coherent rotation because of the inter-slab interaction. The 
SMS saturation field h c s is now higher because its magnetization is stabilized by the interaction 
with the HMS. On the other hand, h c h = h sa .t decreases because the molecular field of the already 
reversed SMS acts against the HMS anisotropy field. For stronger EI (Aei = 0.11,1) the plateau 
disappears completely {h c s = hi = h sat ), indicating that once the SMS magnetic moments reach 



24 



Ah Ah 




-0,06 -0.04 -0,02 0,02 0,04 0,06 

h 

(c) 



Figure 15: (a) Effect of EI, (b) of DDI, and (c) of DMI (along the z axis) at the HMS-SMS interface on 
the hysteresis cycle and switching mechanisms of the HMS/SMS-SI multi-layer systems with Afh = 12 and 
Af s = 16. 

a certain deviation, the EI induces a cascade effect in the HMS that causes the switching of the 
latter. The same observations of nonuniform magnetization switching leading to a smaller coercive 
field were made for magnetic recording exchange-spring media |16[ [T7] . Finally, an increase of the 
EI increases the remanent magnetization. The coercive field increases as long as h sa ,t > h c s . As 
soon as the EI becomes strong enough to observe the cascade effect (/i sa t = h^), the coercive field 
becomes a decreasing function of EI. Now, since the inter-layer exchange coupling decreases for an 
increasing spacer thickness, the coercive field passes by a minimum when the latter increases, again 
in agreement with what has been observed in Refs. [TBJ [TTJ 
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Weak (Addi = 0,0.03) or medium (Addi = 0.1,0.3) DDI has a similar behavior to that of the 



EI, as can be seen in Fig. 15 (b). There is some difference, however, between the EI and DDI 
regarding the evolution of the hysteresis cycle as the inter-slab coupling is varied. While increasing 
EI decreases h sa ,t, the DDI induces an increase of h sa ,t with increasing Addi- Indeed, DDI whose 
bond is along the z axis induces an additional anisotropy at the interface along this axis and thereby 
tends to stabilize the magnetization in this direction. For Addi = 1-0 the induced anisotropy field is 
stronger than the HMS anisotropy field. This implies that it is possible for the magnetic moments 
of the HMS that are far from the interface to switch before those at the interface. When the applied 
field becomes in excess of the induced-anisotropy field a complete switching is achieved. Finally, 
both the remanent magnetization and coercive field increase with DDI. 



The effect of the DMI along the z axis on the hysteresis cycle is shown in Fig. 15 (c). A similar 
behavior to those exhibited by EI and DDI for weak interaction (Admi = 0.09) is observed. Again, 
we observe a progressive rotation of the SMS magnetization before the applied field becomes strong 
enough to induce a reversal of the HMS magnetization. As we discussed earlier, with no applied 
field the DMI along the z axis induces a variation in the deviation of the magnetic moments of the 
HMS, while pinning the SMS (net) magnetic moment along the x axis. The slight deviation induced 
in the magnetic moments of the HMS leads to a reduced HMS switching field h c h = h sa ,f Further 
increase of the interaction (Admi = 0.3, 1) induces again a cascade effect in the HMS that causes 
both slabs to switch in the same field {h c s = hf). However, we have to note that unlike EI and 
DDI, the HMS/SMS switching field is not equal to h sa t- This is again due to the deviation induced 
in the HMS magnetic moments by DMI. Both the remanent magnetization and the coercive field 
decrease with increasing DMI. 

We can see that each kind of interaction induces a specific "single magnetic moment"-like behavior 
for strong coupling A. A strong EI tends towards a system with no anisotropy, whereby the area of 
the loop tends to vanish, reaching saturation with very low applied fields in a switch-like behavior. 
DDI on the other hand tends to induce a square loop, typical of a very high uniaxial anisotropy 
with easy axis parallel to the applied field. Finally, the DMI tends to narrow the cycle and the 
magnetization follows the applied field, typical of a system with uniaxial anisotropy perpendicular 
to the field. Obviously, the perfect "single magnetic moment" behavior is never reached because 
the inter-slab interaction is present only at the interface and the deviation of the outer sub-layers 
is limited only by the intra-slab exchange interaction and the applied field. 
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V. CONCLUSION 

We have studied a magnetic multi-layer system composed of a hard magnetic slab with out-of- 
plane uniaxial anisotropy and a soft magnetic slab with in-plane uniaxial anisotropy, separated by 
a nonmagnetic spacer. We have considered three cases of inter-slab coupling, namely exchange, 
dipolar or Dzyaloshinski-Moriya interactions. The soft magnetic slab has been modeled as a stack 
of atomic (e.g. Fe) layers, while the hard magnetic slab has been modeled either as a macroscopic 
magnetic moment or as another stack of atomic (e.g. FePt) layers. Each atomic layer is modeled as 
a macroscopic magnetic moment representing its net magnetic moment, and is coupled to adjacent 
layers by exchange coupling. 

We have investigated the effect of the external magnetic field, the in-plane or out-of-plane 
anisotropy, and the three interactions on the deviation angle (relative to the hard slab anisotropy 
easy axis). We have computed the magnetization profile through the whole multi-layer system. 
Our computing method consists in solving the set of (coupled) Landau-Lifshitz equations for the 
net magnetic moments of the layers. We first validate this method by comparing the corresponding 
results to the previously obtained analytical expressions for the case of exchange inter-slab coupling. 

For the effect of the applied magnetic field, we found that with sufficient number of layers in 
the soft slab, the system behaves according to the Stoner-Wohlfarth regime and that there exists 
a critical field at which the whole system aligns along the applied field. 

For exchange and dipolar interactions, there is an asymptotic value that depends on the 
anisotropy and the intra-layer exchange. For the dipolar interaction with a bond along the hard 
slab anisotropy, this asymptotic value is given by the analytical expression for rigid interface, where 
the hard slab is modeled as a single pinned magnetic moment, and only the variation in the soft 
slab is relevant. The dipolar interaction was next extended through the whole soft slab with a rigid 
interface. The ensuing effect on the magnetization profile of the soft slab has been recovered by 
an effective dipolar coupling at the interface only, upon re-scaling the intra-slab exchange coupling. 
The two corresponding effective coupling parameters depend on the initial dipolar interaction, but 
not on the number of SMS layers. 

For the Dzyaloshinski-Moriya interaction, we found that the magnetic moments of one of the 
slabs are pinned in a given direction, whereas those of the remaining slab rotate in either the polar 
angle or the azimuthal angle, depending on the direction of the vector D. Large values of the 
DM coupling lead to a system with rigid interface, with the magnetic moments at the interface 
being perpendicular to each other and to the vector D. In addition, a switch-like mechanism can 
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be achieved with this interaction. Indeed, an indirect reversal of the SMS magnetization can be 
achieved by directly forcing a reversal of the HMS magnetization with the help of a magnetic field. 
It is obviously also possible to obtain the desired effect for the exchange-spring system by achieving 
the reversal of HMS via the switching of the SMS magnetization with the help of a smaller magnetic 
field. 

A comparison between the three interactions with typical orders of magnitude has been given. 
The exchange coupling shows the strongest effect, and when added, the dipolar and Dzyaloshinski- 
Moriya interactions induce a slight (but non negligible) deviation in either the polar or azimuthal 
magnetization profile. 

Finally, hysteresis cycles for the three different interactions are computed. A typical exchange 
spring behavior, where the soft slab switches first, followed by the hard slab at a stronger field, is 
observed for weak interaction in all cases. Strong coupling causes both slabs to switch under the 
same field. Exchange and Dzyaloshinski-Moriya interactions tend to narrow the cycle, while the 
dipolar interaction leads to squared cycles. For application purposes such as magnetic recording 
using vertical exchange-spring media, the Dzyaloshinski-Moriya inter-layer coupling strongly re- 
duces the coercive field while keeping high values of the anisotropy which thus ensure good thermal 
stability. 

A work in progress consists in treating each atomic layer as a two-dimensional lattice with the 
aim to compute spin correlations as functions of the various energy parameters and to determine 
the spin-wave spectrum. In the near future, this zero-temperature study will be extended to finite 
temperature with the aim to investigate thermal effects with special emphasis on the calculation of 
activation rates of such multi-layer systems and thereby assess their thermal stability. 
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